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• • . 1 INTRODUCTION 



ABSTRACT 

Using high resolution (~85,000) and high signal-to-noise ratio (~200) optical spec- 
tra acquired with ESO/UVES spectrograph, we determined interstellar column den- 
sities of C2 for six Galactic lines of sight with E(B-V)'s ranging from 0.33 to 1.03. 
For the purpose we identified and measured absorption lines belonging to the (1,0), 
(2,0) and (3,0) Phillips bands A 1 !!, - X 1 !^. Identification of a few lines of the C 2 
(4,0) Phillips system towards HD 147889 is reported. The curve of growth method 
is applied to equivalent widths for the determination of column densities of individ- 
ual rotational levels of C^. Excitation temperature is extracted from the rotational 
diagrams. Physical parameters of the intervening molecular clouds: gas kinetic tem- 
peratures and densities of collision pa rtners were estimated from c omparison with the 
theoretical model of excitation of Ci (jvan Dis hoeck fc B lack 1983). 
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Diatomic carbon is found in a va riety of celestial bodie s 
such as comets (C2 Sw an band - iMaver fc O'delll (|l968l ); 
ILambert fc Danka Jl98&f); Ballik- Ramsay a nd Phillips bands 



I Johnson et al 



"(ll983l);lGredel et all { 19891)) the S un (Swan 



band - Greves se fc Sauva ll 1 19751 ) ; iLamberd (|l978l )) 



bon stell ar atmospher e s (Sw an, Ballik-Ramsay and Phillips 
bands - iQuerci et all (|l97ll )), circumstellar shells of car- 
bon rich po s t-AG B stars (Phill i ps an d Swan bands - 
iBakker et ail |l99rj ); uBakker et all |l997l )) and in interstel- 
lar clouds. 

In the interstel lar medium, Ci was first detected by 
ISouza fc LutJ (|l977h . They observed absorption lines of the 
(1,0) Phillips band A 1 II U -X 1 E+. Lines of the (2,0) Phillips 
sy stem were first tentati vely detected in interstellar clouds 
bv lChaffee fc Lutzl (|l978h . The evident demonstration of the 
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band was shown by IChaffee et all dl98Ch. The Phillips ban d 
(3,0) was discovered only by van Dishoeck fc Black! (l986). 
Some C% lines of the (4, 0) Phillips b a nd in direction to 
HD 204827 were shown bv lHobbs et~afl (|2008l ) (few of very 
weak and difficult to measure lines of this band are demon- 
strated here for the first time toward HD 147889). 

C2 is particularly interesting because it is the simplest 
multi-carbon molecule. Its abundances give information on 
the chemistry of interstellar clouds, especially on the path- 
way to the formation of long chain carbon molecules which 
may be connected with carriers of diffuse int erstellar bands 
(DIBs) (|Douglaslll977l : iThorburn et al.ll2005t ). Additionally, 
the analysis of C2 lines allows to determine physical condi- 
tions in interstellar clouds. 

In the interstellar medium, because of the very rare 
collisions and low temperatures, we expect only absorption 
lines from the groun d elect r onic s tate. For Ci, there are Mul- 
liken (discove red bvlSnowl (Il978l). D^ Y± - X l E+, ~2313A, 



see also - ILambert et al 



Sonnentrucker et al.l 



l|2007h ). F-X (first detected bv lLienl l| 19841 ). F 1 !!,, 



6900- 



T342A) and Phillips systems (A 1 n il 
12000 A). 

C2, as a homonuclear diatomic molecule, has a neg- 
ligible dipole moment and hence radiative cooling of the 
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excited rotational level s may so only through the slow 
quadmpole transitions (|van Dishoeck fc Black! Il982h . The 
rotational levels are pumped by the galactic interstellar ra- 
diation field and excited effectively above the gas kinetic 
temperature. The rotational ladder of the electronic absorp- 
tions from the high rotational levels (here up to J" =26) 
are usually observed. Because of that, lines of the diatomic 
carbon from a long-lived ground state rotational levels are 
measurable and can be the sensitive diagnostic probes of 
conditions in molecular clouds that produce the interstellar 
absorption lines, in contrast to polar molecules, such as CH 
or CN, where usually only a few absorption lines from the 
lowest rotational levels are observed. 

Relative abundances of C2 were predicted for the in- 
terstellar clouds on the basis of detailed chemical models 
(|Black fc Dalgarno|[l977l ; iBlack et alJll978l ). The excitation 
mechanisms of C2 have been already analysed in detail 
|Chaffee et al.lll980l ; lvan Dishoeck fc Blackll 19821 ). 

The main purpose of this paper is to determine basic 
physical parameters like temperatures and densities of the 
intervening clouds. Four out of the six analysed objects (Ta- 
ble 1) were stud ied before, but high quality spectra from 
the ESO archive (|Bagnulo et al. 2003) allowed to accurately 
measure these weak features. Broad spectral range of the 
UVES spectra allows us to analyse up to four bands by con- 
trast with previous analyses usually based on single (2,0) 
band of the Phillips system. Simultaneous observations of 
different bands give possibility to compare individually de- 
termined results for each band and make results more reli- 
able. In this paper two new objects with the interstellar ab- 
sorption lines of C2 are presented, increasing the presently 
av ailable sample of interstella r clouds (24 lines of sight - 
see lSonnentrucker et alj |2007l ') - Table 13) where a detailed 
analysis of excitation of C2 was made (estimates of densities 
and excitation temperatures). 

The next section describes our data, the observations, 
reduction and criteria for choosing stars. In Sect. 3 we intro- 
duce methods of analysis of the observational data. General 
discussion and summary of our conclusions are given in Sect. 
4 and 5. Discussion of individual results for each star of the 
sample and comparison between our results and those of 
previous papers is demonstrated in the Appendix. 



2 THE OBSERVATIONAL DATA 

We used the archived spectrosc opic observations co llected 
using the spectrograph UVES (|Dekker et al.l l2000h . UV- 
Visual Echelle Spectrograph (UVES) is the high-resolution 
spectrograph of the VLT fed by the Kueyen telescope of the 
ESO Paranal Observatory, Chile. UVES is a cross-dispersed 
echelle spectrograph designed to operate with high efficiency 
from the atmospheric cut-off at 3000 A to the long wave- 
length limit of the CCD detectors (about 11,000 A) so it 
is really suitable instrument allowing to observe four bands 
of the Phillips system in single exposure. The high-quality 
of UVES allows to get excellent spectra with the maximal 
resolution of about 85,000 (accurate value in Table 1) and 
signal-to-noise ratio above 200. 

We have chosen spectra with interstellar molecular lines 
of C2 towards six early-type stars (Table 1) from the UVES 
archive. Objects with intermediate colour excess and with 



Table 1. Basic data for program stars (Thorburn et 
al. 2003, Hunter et al. 2006). 



object 


name 


Sp/L 


V 


E(B-V) 


R 


HD 76341 




091b 


7.17 


0.49 


79,000 


HD 147889 




B2V 


7.90 


1.07 


98,000 


HD 148184 


X Oph 


B2V 


4.28 


0.52 


81,000 


HD 163800 




04V 


7.01 


0.60 


82,000 


HD 169454 


V430 Set 


BHa 


6.65 


0.93 


110,000 


HD 179406 


20 Aql 


B3V 


5.34 


0.33 


63,000 



one really dominant velocity component, at the resolution 
of observational material, were selected to make it likely that 
C2 molecular lines originate from single clouds. In addition 
KI 7699 A line profiles were checked for possible existence of 
more than one dominating Doppler components. Generally 
we cannot exclude exis tence of multiple closely - spaced com- 
ponents. For example, IWeltv fc Hobbsl (|2001 ) show a very 
weak Doppler component in KI and CH toward HD 148184 
which are not visible in our spectra. Presence of weak com- 
ponents should not contaminate rather weak C2 lines. These 
authors also show that almost all of the program objects 
have multiple components in Nal; ho wever this componen t 

~Jl2007l) 



may be not related to molecular ones IjBondar et al 



In ult ra-hi gh resolution spect r a of HD 169454, Crawford! 



l|l997l ) and lCrawford fc Barlowl (|l996l ) see two Doppler com- 
ponents in C2 lines. These are yet unresolved at the UVES 
resolution. 

The spectral an alysis was made with the Dech20T code 
|Gafazutdinovlll992f ). 

Spectral regions with interstellar C2 absorption lines, 
especially the (3,0) and (4,0) bands, are contaminated by 
atmospheric lines. Our spectra were divided by spectrum of 
proper standard (Spica) to remove telluric lines. In this way 
a majority of telluric features was removed. Remaining lines 
were cut out one by one manually. The continuum was traced 
to remove broad stellar absorption features. Also possible 
broad DIB at 7721.7 A (|Herbig fc Lekalll99ll ) blended with 
C2 (3,0) Q(2) line was removed in this procedure. 

The final spectra were normalised to unity to en- 
able measurements of the equivalent widths. The equiva- 
lent widths were measured by fitting Gaussian profile to 
each absorption line. At the resolution of the spectra (~ 
3 — 5kms~ 1 ) the profiles of single molecular lines would 
have the shape of the instrumental profile. Expected Doppler 
broadening of the profiles c aused by therma l and turbulent 
motions (< 0.7kms _1 - see lCrawfordl (Il997t ) (Table 4.)) are 
significantly lower. 

The uncertainties of the equivalent widths (AEW) were 
estimated with the /iMifj] splot task taking into consider- 
ation signal-to-noise ratio in the portion of spectrum close 
to measured C2 line. In these cases, where continuum trac- 
ing was uncertain, additional uncertainty was estimated by 
changing level of continuum. The errors were propagated to 



The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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Figure 1. The regions of the Ci (1,0), (2,0), (3,0) and (4,0) Phillips band in the spectrum of HD 147889. Spectrum was normalised to 
a continuum level of 1. C2 absorption lines are indicated. The wavelength scale has been shifted to the rest wavelength velocity frame 
using the interstellar KI line (7698.965 A). 
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the uncertainties of the determined parameters (e.g. excita- 
tion temperatures, column densities) and used in search of 
the best fit model parameters. 

In summary, we measured absorption lines (P, 
Q, R branches) in bands (1,0) 10133 - 10262 A, (2,0) 
8750 - 8849 A, (3,0) 7714 - 7793 A. In one case, HD 147889, 
we were able to identify and measure interstellar absorption 
lines of the (4,0) 6909 - 6974 A band (Fig. 1). 

The final spectra of HD 147889, with the continuum 
level normalised to unity, are shown in Figure 1. Broad 
stellar lines were removed in the procedure of continuum 
tracing. There are plenty of weak absorption lines (Fig. 1) 
thanks to the high quality of the UVES data. Even the P(4) 
and Q(8) lines, which usually are blended in the (2,0) band, 
are well separated in all objects. A lot of C2 lines were also 
found in the (3,0) band. The oscillator strengths for (4,0) 
are about five times weaker than those in the (2,0) and, in 
addition, molecular lines around 6915 A are seriously con- 
taminated by telluric lines. C2 lines of the Phillips system, 
seen in all spectra of the program stars, allow to identify 
rotational components of P, Q and R branches to high ro- 
tational levels, up to J" = 16. In the spectrum of HD 147889 
we managed to identify Q components up to J" =26. 



3 RESULTS AND INTERPRETATION 

We have identified lines of the (1,0), (2,0), (3,0) and in one 
case (4,0) bands of the C 2 Phillips system (A x n u - X X E+). 
The equivalent widths with errors of all measured interstel- 
lar lines of C'2 towards six program stars are given in Ta- 
bles 2-5. Poor quality measurements of equivalent widths 
are marked with sign u behind the value. These uncertain 
values were not used in the further analysis of the column 
densities. 

Column density of a rotational level J" may be derived 
from the equivalent width W^mAlofthe single absorption 
line using the relationship (Frisch 19721 ) 



N col = 1.13 x 10 1 



(1) 



where A is the wavelength in [A] , /y is the absorption oscil- 
lator strength. 

The energies of the lower rotati onal level were de- 
termined using molecular constants of M arenin fe Johnson! 
(1970). The wavelengths are generally dete r mined from 
laboratory wavenumber s of IChauville et al.l (|l977T l and 
iBallik fc Ramsavl l| 19631 ) con verted to air w avelengths us- 
ing Edlen's formula following iMortonl l|l99ll ). Wavelengths 
of three lines R(2) P(2) and P(4) of the (2,0) band, absent 
111 (IChauville et al.ll 19771) . were computed with 
19881 ) spectroscopic constants. According to 



Douav et al 



Douav et al 



the line positions calculated with their constants 
should be more accurate than the previous measurements. 
The oscillator strengths correspond to vibrational oscilla- 
tor strengths f 10 = 2.38 x 10" 3 , /20 = 1.44 x 10" 3 , /30 = 
6.67xl0" 4 , /jo = 2.71 xlO" 4 . The oscillator strengths for in- 
di vidual transiti o ns we re computed according to description 
in iBakker et all j 19961 ) using their code M OLLEY. Vibra- 
tional oscillator strengths we re taken from Langhoff et al.l 
|l990l ) for (1,0) and (2, 0). IBakker et all ll 19971) (citing 
Langhoff) (3,0) and from Ivan Dishoeckl (| 1983) for (4,0). 



Source of band origins were IChauville et~afl l|l977l ) and 
IBallik fc Ramsavl j 19631 ). 

Equation (1) is accurate for the optically thin case 
(when the absorption lines are on the linear part of the 
curve of growth). Only a few lines from the (1,0) band 
of HD 147889, where the equivalent widths are the largest, 
are evidently optically thick. Then curve of growth method 
was applied for the derivation of column densities instead 
of the equation (1). For an optically thick line we deter- 
mined the turbulent velocity through the minimalization 
of the dispersion of column densities for each level. We 
checked various values of the velocity dispersion parame- 
ter (b = 0; 0.5; 1; 1.5 kms -1 ) and 0.5 kms -1 was found to 
give the lowest dispersion. This value was then applied to all 
of the program stars. It is consistent with the value derived 
in oth e r studies o f mole cular absorptions (e.g. iGredel et all 
|l99ll) ; ICrawfordl (|l997f )V The application of equation (1) 
to optically thick lines underestimates column densities by 
28 per cent in the worst case of Q(2) (1,0) absorption line 
in HD 147889. 

The resulting column densities for each rotational level 
N co i(J"), derived uncertainties and the number of measure- 
ments used to determine that value are presented in Table 6. 
The total C2 column density, defined as the sum of the mean 
column densities of the observed levels and of the contribu- 
tion of the unobserved levels estimated from the theoretical 
model characterised by the best-fit parameters (see below), 
are shown in column 5 of Table 7. There were available three 
Phillips bands so the column density could be maximally es- 
timated from 9 measurements. Not every part of spectrum 
allows to measure very weak interstellar features of C%. In 
some cases column density was derived from only one tran- 
sition (e.g. for the highest rotational levels J" = 10, 12, 14, 
16) . Absorption lines from lower levels (J" = 2 and 4) are the 
most populate d and the easiest to be measured accurately. 

Following Ivan Dishoeck fc Black! (|1982| ) we present 
derived column densities in form of rotational dia- 
grams (Fig. 2) where weighted relative column densities 
-ln[5N co i(J")/(2J" + l)A co; (2)] are plotted versus energy 
of lower level E" /k (where E" is the energy of the rotational 
level J" and k is the Boltzmann constant). The errorbars cor- 
respond to the derived uncertainties. The slope of a straight 
line on this diagram is nicely connected to the excitation 
temperatur e, a = —1/T exc . It is well kn own from previous 
works (e.g. Ivan Dishoeck fc Black! ||1982|) ) that populations 
of all rotational levels cannot be characterised by a single 
rotational temperature. The lowest J" levels are described 
by the lower excitation temperature than higher levels. Such 
behaviour of the rotational le vels was nicely described in th e 
model of excitation of C2 by Ivan Dishoeck fc Black! l|l982l ). 
In this model, the molecule is heated by the electronic tran- 
sitions from the ground state, and subsequently cooled down 
through cascading to the ground electronic state through ex- 
cited vibrational levels. Because C2 is a homonuclear species, 
the mechanism of cooling is inefficient and the high rota- 
tional levels are significantly populated. The population of 
the lowest rotational levels is influenced by the collisions 
with the gas, mainly atomic and molecular hydrogen (hence 
the density of collision partners n c — uh + nH 2 )- 

For the interpretation of the rotational diagrams we 
have constructed a grid of models based o n the radiative 
excitation model of Ivan Dishoeck fc Black! l| 19821 ). The de- 
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Table 2. Observation summary table with equivalent widths [mA] of Ci (1,0) Phillips lines toward program stars. In table are also 
B(N" = J") - branch identification (J" - low rotational level) and A - wavelength in air in A (see text for the references). 



B(N" = J") 


HD76341 


HD 147889 


HD148184 


HD163800 


HD169454 


HD179406 


A [A] 


P(6) 




8.3 ± 1.5 


1.6 ± 1.0 


0.6 ± 0.5m 


4.0 ± 1.0 


3.2 ± 1.0 


10133.603 


R(8) 




5.8 ± 1.5 


2.4 ± 1.5 


1.0 ± 0.7 


0.7 ± 1.3 


1.9 ± 0.8m 


10133.854 


Ti( A \ 

R{4) 


2.2 ±1.1 


15.6 ± 1.4 


6.5 ± 1.6 


3.6 ± 1.3 


7.0 ± 1.2 


8.7 ±1.1 


i A-i or 1 a r\ 

10135.149 


R(10) 




3.4 ± 2.2 




0.5 ± 0.5m 




2.0 ± 1.4m 


i a-i or / \ o o 

10135.923 


R(2) 


3.2 ± 1.2 


17.1 ± 1.4 




2.6 ± 1.0 


12.5 ± 1.1 


6.1 ± 1.5 


10138.540 


R(12) 






0.5 ± 0.9 








1 n 1 on or\tr 

10139.805 


R(0) 


1.6 ±1.0 


12.4 ± 1.4 


4.0 ± 1.2 


3.8 ± 1.0 


14.3 ±1.0 


5.2 ±1.0 


i r\i A O TOO 

10143.723 


«(14) 








0.8 ± 0.5 






1U145. DUD 


Q{2) 




22.1 ± 1.4 


6.4 ± 0.5 


4.3 ± 1.3m 


21.0 ± 1.4 


9.0 ±1.3 


10148.351 


0(4) 


2.6 ± 1.2 


on q _l_ i a 
2U.3 ± 1.0 


7.2 ± 1.2 


3.2 ± 1.4m 


12.5 ± 1.0 


11.4 ± 1.2 


lU10l.02o 


P{2) 




5.6 ± 1.4 




2.0 ±0.9 


3.9 ± 1.2 




10154.897 


0(6) 




17.1 ± 1.4 


4.4 ± 1.3 


6.0 ± 1.4m 


5.9 ± 1.0 


6.3 ±1.0 


10156.515 


0(8) 


0.7 ±0.7 


13.7 ± 1.2 


2.8 ± 1.5 


4.6 ± 1.3 


4.5 ± 1.2 


4.2 ±1.2 


10163.323 


P(4) 


0.3 ± 0.3m 


8.8 ± 1.5 


2.1 ±0.9 




2.8 ± 1.0 




10164.763 


Q(10) 


0.8 ±0.8 


10.6 ± 1.6 


2.6 ± 1.2 




1.7 ±0.8 


3.0 ± 1.0 


10171.963 


P(6) 




8.8 ± 1.5 


1.2 ±0.9 




3.0 ±0.9 




10176.252 


0(12) 




4.7 ± 1.2 






0.9 ± 0.9m 




10182.434 


P(8) 




4.9 ± 1.3 


1.4 ±0.7 




1.3 ±0.6 




10189.693 


0(14) 




3.3 ± 1.2 










10194.755 


P(10) 




2.8 ± 1.6 










10204.998 


Q(16) 




4.0 ± 1.6 










10208.931 


P(12) 




1.8 ±0.9 










10222.171 


P(14) 




1.7 ± 0.8 










10241.246 


Q(18) 




1.8 ± 1.4m 










10224.982 



Table 3. The same as in Table 2. but for the (2,0) Phillips band. Wavelengths marked with 1 were computed from spectroscopical 
constants of Douay ct al. 1988. 



B{N" = J") 


HD76341 


HD 147889 


HD148184 


HD163800 


HD169454 


HD179406 


A[A] 


R(6) 


0.6 ±0.4 


6.1 ±0.3 


0.7±0.4 


1.6 ±0.5 


2.4 ±0.3 


1.0 ±0.5 


8750.847 


R(8) 




3.6 ±0.3 


0.6 ±0.5 


1.0 ±0.4 


0.7 ±0.3 


0.4 ±0.4 


8751.487 


fl(4) 


0.9 ±0.5 


9.2 ±0.3 


1.5 ±0.4 


2.0 ±0.5 


3.7 ±0.3 


2.7 ±0.5 


8751.684 


R(10) 




2.2 ±0.3 


0.7±0.4 




0.4 ± 0.3m 




8753.578 


R(2) 


1.4 ±0.5 


10.5 ± 0.3 


2.7±0.4 


1.8 ±0.4 


8.1 ±0.3 


4.1 ±0.5 


8753.945 1 


R(12) 




1.8 ±0.3 






0.3 ±0.2 




8757.127 


R(0) 


1.1 ± 0.4 


7.5 ±0.3 


1.4 ± 0.4m 


1.8 ±0.4 


8.5 ±0.3 


3.3 ±0.4 


8757.683 


Q(2) 


1.2 ± 0.5 


12.7 ±0.3 


3.0 ±0.4 


2.7 ±0.5 


10.3 ±0.3 


4.7 ±0.4 


8761.194 


R(14) 












0.9 ± 0.3m 


8762.144 


Q(4) 




13.1 ± 0.3 


3.3 ±0.5 


3.5 ±0.7 


7.9 ±0.3 


7.1 ± 0.5m 


8763.751 


P(2) 




2.4 ±0.3 


0.5 ± 0.5m 


0.5 ±0.4 


2.0 ±0.3 


1.4 ±0.4 


8766.026 1 


Q(6) 




9.2 ±0.3 


3.4 ±0.4 


2.1 ±0.5 


3.0 ±0.3 


3.6 ±0.4 


8767.759 


i?(16) 














8768.628 


Q(8) 


0.3 ± 0.3 


6.3 ±0.3 


1.7±0.4 


1.2 ±0.5 


2.3 ±0.3 


1.3 ± 0.8m 


8773.220 


P(4) 




4.7 ±0.3 


1.4 ±0.5 


0.8 ±0.4 


1.7±0.3 


0.9 ±0.5 


8773.422 1 


Q(10) 


0.5 ± 0.4 


4.5 ±0.3 


1.0 ±0.5 


0.9 ±0.4 


0.8 ±0.3 


2.2 ± 0.5m 


8780.141 


P(6) 


0.6 ± 0.5 


5.1 ±0.3 


1.0 ±0.4 




0.9 ±0.8 


1.7 ± 0.5m 


8782.308 


0(12) 




2.3 ±0.3 


0.6 ±0.4 




1.0 ±0.3 




8788.558 


P(8) 




2.4 ±0.4 






0.8 ±0.3 




8792.649 


0(14) 




2.1 ±0.3 






0.6 ±0.4 


0.6 ± 0.3m 


8798.459 


P(10) 




2.2 ±0.3 










8804.499 


0(16) 




2.2 ±0.3 


0.8 ±0.4 




0.6 ±0.4 




8809.842 


P(12) 




0.8 ±0.4 










8817.827 


0(18) 




1.8 ± 0.5m 










8822.725 


0(20) 




2.2 ± 0.7m 










8837.119 


0(22) 




0.7 ± 0.4m 










8853.041 


0(26) 




0.7 ± 0.5m 










8889.532 
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Table 4. The same as in Table 2. but for the (3,0) Phillips band. 
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/ / lO.OZo 
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771 7 AfiQ 
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771 O QOO 
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10-1-04 
l.U IE U.4 


U.4 El U.41i 








7790 7AS 
/ / ZU. / 4(5 




n 7 -t- o 
u. / in u.ou 


A R 4- A 


9 k 4- 7 
z.o ie u. / 


1 .0 31 u.u 


r: n 4_ n r: 
o.y EE U.O 


174-n Q,, 
1. / EE u.yu 


7799 fiQ^ 

/ / zz.uyo 




n 7 -t- n 

U.I X U.^iU 


A ^ 4- n A 
4.0 IE U.4 


i n 4- n r 

l.U EE U.O 


o 7 4- O Q 
Z. ( EE u.y 


994.(14 
O.Z EE U.4 


o o 4_ n Qi/ 


779A 91 Q 

i i z^.ziy 






U.O EE U.O 










779c: 9/in 
1 ( ZO.Z^iU 


P(2) 










2.0 ±0.8 




7725.819 


0(6) 




4.0 Ei 0.4 






1.4 ±0.6 


2.3 ± 0.5m 


7727.557 


fl(16) 














7730.963 


P(4) 




2.2 Ei 0.6 


0.4 ± 0.4m 




0.8 ±0.5 




7731.663 


0(8) 




2.4 ±0.6 




2.1 ±0.6 




1.1 ± 0.4 


7732.117 


Q(10) 




1.9 ei 0.4 






0.9 ±0.5 




7737.904 


P(6) 




1.8 ±0.4 






0.4±0.4w 




7738.737 


0(12) 










1.0 ±0.5 


0.7 ± 0.7m 


7744.900 


P(8) 






0.3 ±0.3 




0.8±0.5u 




7747.037 


0(14) 














7753.141 


P(10) 




0.5 ±0.4 


0.5 ± 0.5 








7756.582 


0(16) 










0.6±0.6u 




7762.623 


P(12) 














7767.369 



Table 6. C2 column densities N co i(J") [10 12 cm 2 ] for each low rotational level J", (N - number of measurements 
used to determined given N co i). 





HD 76341 




HD 147889 




HD 148184 




HD 163800 




HD 169454 




HD 179406 




J" 




N 




N 


N c „i 


N 




N 




N 




N 





0.9 ±0.3 


2 


8.5 ±0.4 


3 


1.6 ±0.3 


2 


1.9 ±0.3 


2 


9.5 ± 0.4 


2 


3.1 ± 0.4 


2 


2 


3.4 ± 0.6 


3 


30.7 ±0.6 


8 


6.7 ±0.4 


5 


4.9 ±0.6 


6 


23.8 ±0.6 


6 


10.4 ±0.7 


6 


i 


2.8 ±0.7 


3 


31.3 ±0.6 


9 


6.8 ±0.6 


8 


6.4 ±0.9 


6 


14.8 ±0.5 


6 


10.8 ± 1.0 


1 


6 


2.0 ± 1.4 


1 


22.1 ±0.6 


9 


5.0 ±0.6 


7 


4.7 ±0.9 


2 


6.8 ± 0.5 


6 


6.4 ±0.6 


4 


8 


0.6 ±0.5 


2 


13.7 ±0.5 


8 


3.2 ±0.6 


6 


2.7 ±0.7 


4 


4.0 ±0.5 


6 


3.4 ±1.1 


2 


10 


0.9 ±0.6 


2 


9.6 ±0.5 


9 


2.3 ±0.7 


1 


1.9 ±0.8 


1 


1.8 ± 0.5 


3 


2.9 ± 1.0 


1 


12 






5.2 ±0.5 


6 


1.1 ±0.7 


2 






1.8 ±0.5 


3 






14 






4.1 ±0.6 


1 






1.3 ±0.9 


1 


1.2 ± 0.8 


1 






16 






4.7 ±0.6 


2 










1.2 ± 0.8 


1 







Table 5. The same as in Table 2. but for the (4,0) Phillips band. 



B(N" = J") 


HD147889 


A[A] 


i?(4) a 


0.5 ± 0.5m 


6909.412 


R(2) 


3.3 ± 1.4m 


6910.577 


0(2) 


2.0 ± 1.0m 


6914.975 


0(4) 


1.2 ± 0.8m 


6916.788 


0(6) 


1.4 ± 0.8m 


6919.656 



a R(4) may form a blend with R(6) 6909.374 1 



tailed behaviour of excitation temperature depends on the 
gas kinetic temperature and on the ratio of the collisional 
rate n c x a to intensity of the average galactic field, here 
expressed by I. The parame ter I has been introduced by 
Ivan Dishoeck fc Black! |l982) as a scaling factor of the stan- 
dard field adopted in their paper. With the lack of detailed 
information on the radiation field in individual diffuse clouds 



we made a crude assumption 1 = 1. Then, from the best fit- 
ting models of excitation, one can estimate the density of 
collision partners in a molecular cloud n c . 



Following Ivan Dishoeck fc Black! ( |l982 ) , we computed 
grid of models with step of 5 K in kinetic temperature and 
25 cm -3 in collisional partner densi ties. The models were 
constr ucted by solving eqn. 27 in Ivan Dishoeck fc Black 
(1982), using their quadrupole and collisional rates and 
radiative excitation matrix. To adjust for the differ- 
ences in adopted f-values of Phillips transitions we 
rescaled the density of collisio nal partners according to 
Ivan Dishoeck fc de Zeeuw (1 19841) prescription by fact or 1.59 
(factor of 1.35 from Ivan Dish oeck fc de Zeeuw! l| 19841) . times 
1.18 for difference in f-values used here). The m odels 
with original f-values of Ivan Dishoeck fc Black! (|1982| ) were 
checked to reproduce results of original paper with reliable 
accuracy of 20 percent in the worst case of the higher levels. 
The grid of models was used to find the best fit to the ear- 
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Figure 2. Relative C'2 rotational population diagrams toward six program stars, as a function of the excitation energy (or rotational 
quantum number J"). The solid lines represent fit to the theoretical model, based on the analysis of van Dishoeck & Black (1982). The 
straight dashed line shows the best-fitting Tq2- 



Table 7. Summary of the observational data for C2, T02 (T04, Toe) - rotational temperature calculated from the 
two (three, four) lowest rotational levels, N co i - total column densities and the results of a model: Tki„ - gas kinetic 
temperature, n c - the effective density of collision partners and N co i(J" = 2) - column density derived from J" = 2. 



object 


T 2[K] 


T 0i [K\ 


Toa[K] 


AT coi [10 12 cm- 


" 2 ] T kin [K] 


n c [cm 3 ] 


JV co! (J" = 2)[10 12 cm- 2 ] 


HD 76341 


48 ±48 


47 ± 16 


52 ± 15 


11 ± 1 


34 ±20 


300 ± 200 


3.3 ±0.5 


HD 147889 


49 ±7 


62 ±3 


71 ±2 


133 ± 1 


39 ± 2 


199 ± 7 


30.8 ±0.5 


HD 148184 


82 ± 82 


65 ± 12 


74 ±7 


30 ±3 


45 ± 12 


209 ± 60 


6.6 ±0.4 


HD 163800 


24 ±8 


66 ± 17 


76 ± 13 


28 ± 1 


38 ± 15 


169 ± 43 


5.6 ±0.5 


HD 169454 


23 ± 2 


31 ± 1 


36 ± 1 


65 ± 1 


19 ± 2 


326 ± 17 


23.8 ±0.5 


HD 179406 


38 ± 12 


59 ±9 


62 ± 5 


39 ± 1 


38 ± 9 


251 ± 83 


10.8 ±0.5 



lier determined column densities individually for each tran- 
sitions weighted by their errors. We decided to search for the 
best fit to the absolute column densities, instead of the rel- 
ative to J" = 2 populations. Hence, one additional parame- 
ter of the model was absolute column density N(J" — 2) of 
J" = 2 level (which is equivalent to the total column den- 



sity of C2). This approach is still consistent with the fact 
that kinetic temperature and density of collisional partners 
depends only on the relative populations but changes the 
way how observational errors are propagated in the process 
of finding the best fit model. In result, three independent 
parameters: gas kinetic temperature (Thin), collisional part- 
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C2 Phillips (2-0) 




HD 169454 



r _ rrT 



8760 8762 8764 8766 8768 
Wavelength [A] 

Figure 3. The comparison of intensities of few lines of the (2,0) 
band for 3 objects; the spectra were normalised to the line Q(4) 
8763.74 A (2,0) of HD 147889. The excitation temperatures T 06 
of C 2 in HD 148184, HD 147889 and HD 169454 amount 74, 71 
and 36 K, respectively. 



ners density n c , and column density of the J" = 2 level were 
estimated simultaneously for the set of observed column den- 
sities (see Table 7). The best fitted models are presented on 
the rotational diagrams (see Figure 2, the labels describe 
appropriate values of n c ). The analogous procedure applied 
for the relative populations gives very similar results for T^ n 
and ric, well inside determined errors. The difference in ki- 
netic temperatures determined from both methods is less 
than IK except HD 148184 were it is less than 3K. Note 
that because of nonlinearity of the model the errors of best 
fitted parameters may be very asymmetric. Especially in 
these cases when error is a significant fraction of fitted val- 
ues, e.g. HD 76341 and HD 163800, the uncertainty of T^m 
toward lower values will be less than the uncertainty toward 
higher values. In consequence n c could be much higher to- 
ward both objects. 

We also derived a set of rotational temperatures: T02, 
T04, Toe corresponding to the mean excitation temperatures 
derived from a linear fit to logarithm of column densities 
of the first: two, three and four levels, respectively, start- 
ing from J" = respectively. T02 is the best estimator of 
the gas kinetic temperature, T02 Thin, but generally it is 
not well determined, because there is only one available line 
R(0) in each band absorbing from the J" — level. T04 was 
more often published than Toe, but the latter is usually de- 
termined with better precision. On the other side, both T04 
and Toe depend on collisional density and radiation intensity 
in increasing level. Table 7 summarises all these parameters 



2 - 



2 Excitation temperatures determined in this paper are slightly 
different from those announced in the final version of our paper 



iKaimierczak et al] (2009( 1 The reason of differences was rcanal- 
ysis of errors of equivalent widths. 



4 DISCUSSION 

The C2 molecule in interstellar clouds was observed many 
times, but usually measurements were limited to the 
most easily available to ground-based instruments lines of 
the (2,0) Phillips band. Apart from that measurements 
have usually been done with spectro graphs of moderate 
resolution and signal-to-noise r atio ISnowl 1 19781: iHobbsl 
19791 : iDanks fc Lambertl Il985| _ van Dishoeck fc de Zeeuwl 
19841: iGredel fc Muenchlll986l : Ivan" Dishoeck fc Black! Il989l : 
Federman et al.lll994l ~ 

We measured absorption lines of the interstellar C2 to- 
ward six reddened early type stars. The presence of C2 in 
the sight line toward HD 76341 and HD 163800 has not been 
reported before. 

Generally, excitation temperatures T04 in our program 
stars tend to be above 47 K. One evident exception is 
HD 169454 (31 K). This cloud is definitely different from 
the others. There is the lowest gas kinetic temperature. Van 
Dishoeck & Black (1989) suggested that it is a real translu- 
cent cloud, whereas the other clouds are only diffuse inter- 
stellar clouds. 

The behaviour of excitation temperatures may be seen 
directly by inspection of the portions of spectra around the 
Q(J"=4) transition of the (2,0) band (see Fig. 3). Note, that 
the spectrum was renormalized to the strength of Q(4) ab- 
sorption. It is easily seen, that excitation temperature in 
HD 147889 is almost identical to that of HD 148184 and 
higher than the excitation temperature of HD 169454, which 
is confirmed in the quantitative analysis presented above. 

Simultaneous observations of different vibrational 
bands give possibility to compare individually determined 
column densities for each band. For this purpose data 
collected for HD 147889 are especially useful. The ob- 
served column densities determined individually from bands 
(1,0), (2,0) and (3,0) amount (11.2 ± 0.4) x 10 13 cm -3 , 
(13.3 ± 0.2) x 10 13 cm~ 3 , and (13.2 ± 0.8) x 10 13 cm~ 3 re- 
spectively. The observed column density determined form 
all the bands together amounts (13.0 ± 0.2) x 10 13 cm~ 3 . In 
conclusion, the determination of column densities from only 
one of these bands gives reasonable value of the C2 column 
density; but the possibility of analysing lines of many differ- 
ent bands enhances the accuracy of the results. It excludes 
accidental errors (e.g. cosmic rays, telluric lines, instrumen- 
tal errors) which are impossible to remove when only one 
band is available. 

The interstellar absorption lines of C2 toward 
HD 163800 arise in a single velocity component at a radial 
velocity, Vlsr = 7.2 kms -1 with respect to the local stan- 
dard of rest. This radial velocity and position of the star 
suggest an association of the intervening molecular cloud 
wit h a large cloud of cold a tomic hydrogen first surveyed 
by iRiegel fc Crutcherl (|l972l ) (Riegel-Crutcher cloud) as a 
prominent self-absorption feature in the 21 cm line. The 
distance to the cold atomic cloud is constrained by observa- 
tions of opti cal lines in direction to background stars to be 
125 ± 25 pc l|Crutcher fc Lienlll984al ). Apparent association 
with atomic cloud does not imply physical connection with 
molecular cloud responsible for the observed absorptions of 
C2. Small molecular cloud in the direction to HD 169454 is 
also possibly associ ated with the R iegel-Crutcher HI cloud; 
it was analysed by Ijannuzi et al.l (1 19881 ). The distance to 



HD 163800 from the spectroscopic data is estimated to be 
about 1.5 kpc. 

The detailed comparison between our results and these 
of previous papers is shown in Appendix. 



5 SUMMARY 

From the inspection of absorption lines in direction to 
HD 147889, we have found that individual Phillips bands 
(1,0), (2,0) and (3,0) give consistent results, both for column 
densities and excitation temperatures. However not every 
feature is the same easy to measure in all bands; e.g. the 
P(4) and Q(8) transitions are the easiest to separate in the 
(1,0) band and lines of R branch in (3,0) one (R-branch lines 
of the (2,0) band very often are located near to or in the 
wings of the stellar line HI what makes measurements un- 
certain). In this paper measurements of three Phillips bands 
were used to improve reliability of results due to growing 
number of identified lines and their separation. 

The assumption of the identical radiation field for all 
clouds allows to calculate the density of collision partners 
n c ; the found values range from about 170 to 330 cm -3 (100 
to 500 if errors are included). Absolute values may be differ- 
ent because of the peculiar radiation field and approximate 
treatment of the collisional excitation rate of C2 with H2. If 
the rotational population distributions are useful as diagnos- 
tic tools of the cloud densities, it is essential to have better 
understanding of cross sections of rotationally inelastic col- 
lisions of C2 with H2 ■ Determination of n c is less dependent 
on the accurate value of the gas kinetic temperature. The 
visible transitions of C2 are relatively easily observed, (in 
contrast to the far-UV bands of H2 and very weak lines of 
C3), and provide an important tool for the determination of 
the physical and chemical conditions of diffuse clouds. 

The presence of C2 in the measurable amount is re- 
ported for the two new lines of sight: toward HD 76341 
and HD 163800. The latter star lies in the background of 
the Riegel-Crutcher cloud of the cold HI. The molecular 
cloud, origin of C2 absorptions, may be physically associ- 
ated with the atomic hydrogen cloud, similarly to the earlier 
considered small mo lecular cloud in direction to HD 169454 
|jannuzi et alJll988h . 

To sum up, the high resolution and signal-to-noise ratio 
spectra acquired with the ESO instrument allow us to study 
densities and rotational temperatures varying from object to 
object with increasingly accuracy. This is seen particularly 
well when compare the theoretical model with the observa- 
tional data (e.g. HD 147889 in Figure 2). 

We are planning to continue the survey of high reso- 
lution and high signal-to-noise ratio spectra with detected 
interstellar diatomic carbon and compare with other inter- 
stellar absorption features. It is interesting whether other 
molecules are spatially correlated with C2; this may effi- 
ciently constrain interstellar chemistry. 
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(normalised to continuum, with removed telluric lines) with 
the assignment of the interstellar Ci absorption lines. In 
total we were able to assign 68 absorption features to the P, 
Q, R lines of the Ci Phillips bands. The computed column 
density (13.3 ± 0.1) x 10 13 cm~ 2 is the highest value of the 
sample. 

In the case of HD 148184 we measured 40 absorption 
features belonging to the P, Q, R lines of the C2 Phillips 
bands. 

Toward HD 163800, in total we were able to assign 29 
absorption features to the P, Q, R lines of the C2 Phillips 
bands. We estimated iV coi = (2.8 ± 0.1) x 10 13 cm~ 2 . 

For HD 169494 we were able to identify and measure 50 
absorption lines of the P, Q, R branches of the C2 Phillips 
bands. We estimated N co i = (6.5 ± 0.1) x 10 13 cm~ 2 and n c 
about 330 cm~ 3 (it is the highest value in our sample). 

For HD 179406 we measured 36 absorption features to 
the P, Q, R lines of the C2 Phillips bands. 



APPENDIX A: 

In this section detailed comparison between our results and 
those of previously published is demonstrated. C2 in inter- 
stellar clouds towards HD 147889, HD 148184, HD 169454 
and HD 179406 has already been observed, but towards 
HD 76341 and HD 163800 our determinations are the first 
ones. Comparison between our equivalent widths and data 
from literature is presented in Tables A2-A5. There are some 
differences between our results and those of previous papers 
which could be caused by different quality of the analysed 
spectra. We used the spectra from archival database col- 
lected by UVES, which are characterised by high resolv- 
ing power (A/AA ~ 85, 000) and high signal-to- noise ratio 
(S/N ~ 200). The differences with published data may re- 
semble the difficulty of measuring of weak lines. 

In Table Al there is a comparison between our results 
and those from previous papers. 

The spectrum of HD 76341 was the most noisy in our 
sample, and thus the measurements are less reliable. We 
were able to assign 19 absorption features to the P, Q, R 
lines of the C2 Phillips bands in total. There we found the 
lowest column density (iV co ; = (1.1 ±0.1) x 10 13 cm~ 2 ) of 
C2 in the whole sample. 

In HD 147889 we measured the largest number of C2 
lines. Figure 1 shows portions of the spectrum HD 147889 



C 2 interstellar molecule 



Table Al. The comparison between our results and those from previous papers. Total column densities are rescalcd 
to the value of /20 used in this paper. 



object 


T£- [K] 

kin L J 


T 02 [K] 


T 04 [K] 


T 06 [K] 


Af coi [10 13 cm- 2 ] 


source 




70 


52 ±22 


116 ±28 




8.3 ± 2.1 


Vdll -LVlblUJUCK. 06 Lit; zJticU. w _l yoi 

Sonnentrucker et al. 2007 6 




39 ±2 


49 ± 7 


62 ±3 


71 ±2 


13.3 ±0.1 


this work 


HD 148184 


40 




65 




3.0 ±0.4 
2.4 


van Dishoeck & de Zeeuw 1984 
Danks & Lambert 1983 




50 ± 15 


39 ± 16 


57 ± 12 




2.8 ±0.3 


Sonnentrucker et al. 2007 6 




34 ± 12 


82 ±82 


65 ± 12 


74 ± 7 


3.0 ±0.3 


this work 


HD 169454 


15 

15+10 - 5 
25 








5.8 
4.9 ± 1.0 
9.7 


Gredel & Mucnch 1986 
van Dishoeck & Black 1989 
Gredel 1999 




20 ±5 


24 ±6 


37 ±5 






Sonnentrucker et al. 2007 c 




19 ±2 


23 ±3 


31 ± 1 


36 ± 1 


6.5 ±0.1 


this work 


HD 179406 


55 




57 ±11 


66 d 


3.7 
5.0 ±0.6^ 


Fedcrman & al. 1994 
Sonnentrucker et al. 2007 e 




38 ±9 


38 ± 12 


59 ±9 


62 ± 5 


3.9 ±0.1 


this work 



a gas kinetic temperature is obtained form the best-fitting model. 

b calculated by Sonnentrucker et al (2007) using the equivalent widths of Dishoeck & de Zeeuw(1984). 
c calculated by Sonnentrucker et al (2007) using the equivalent widths of Gredel (1999). 
d excitation temperature from all observed levels Tog 

e calculated by Sonnentrucker et al (2007) using the equivalent widths of Federman et al. (1994). 
f based on equivalent widths of Federman et al. 1994 



M. Kazmierczak et al. 



Table A2. The comparison of our equivalent widths 
[mA] to data of van Dishoeck & de Zeeuw (1984) for 
HD 147889 of the (2,0) band. 



B(N") 


our results 


D&Z a 




6.1 ± 0.3 


4.0 ± 0.5 


n \°) 


3.6 ± 0.3 


^ 2.5 




y.z in u.o 


o. i in u.^± 




z.z in u.o 






l-KJ .O HI U.O 


5.3 ± 0.4 




1.8 it 0.3 






i .o m u.o 


^t.o m u.z 




12.7 ± 0.3 


10.1 ± 0.6 




13.1 it 0.3 


13.9 ± 1.2 




2.4 it 0.3 


2.9 ± 0.3 


Q(6) 


9.2 ±0.3 


7.8 ±0.6 


Q(8) 


6.3 ±0.3 


4.9 ±0.7 


P(4) 


4.7 ±0.3 


4.3 ±0.3 


0(io) 


4.5 ±0.3 


^ 2.5 


P(6) 


5.1 ±0.3 


3.4 ±0.3 


Q(12) 


2.3 ±0.3 




P(8) 


2.4 ±0.4 


^ 2.5 


Q(14) 


2.1 ±0.3 




P(10) 


2.2 ±0.3 




Q(16) 


2.2 ±0.3 




P(12) 


0.8 ±0.4 




Q(18) 


1.8±0.5u 




Q(20) 


2.2±0.7u 




Q(22) 


0.7±0.4u 




Q(26) 


0.7±0.5u 




a Measurements were made with coude echelle spec- 


trograph (with Rcticon) which was fed by the 1.4-m 


coude auxiliary telescope at ESO, La Silla, Chile. 


The resolving power in those observations was 80,000. 


UVES allows to detect more 


transitions due to higher 


resolution and S/N ratio. 




Table A4. The comparison 


of our equivalent widths 


[mA] to data of Fcderman et al. (1994) for HD 179406 


of the (2,0) band. 




B(N") 


our results 


Fedcrman a 


P(6) 


1.0 ±0.5 


3.0 ±0.3 


R(8) 


0.4 ±0.4 




fl(4) 


2.7 ±0.5 


3.2 ±0.2 


R{2) 


4.1 ± 0.5 


5.0 ±0.4 


R(0) 


3.3 ± 0.4 


2.0 ±0.2 


0(2) 


4.7 ±0.4 


7.3 ±0.5 


i?(14) 


0.9±0.3u 




Q(4) 


7.1 ± 0.5u 


5.7 ±0.5 


P(2) 


1.4 ±0.4 


sg 1.1 


Q(6) 


3.6 ± 0.4 


4.2 ±0.4 


Q(8) 


1.3±0.8u 


1.1±0.3 6 


P(4) 


0.9 ± 0.5 


1.8±0.3 f ' 


Q(10) 


2.2±0.5u 




P(6) 


1.7 ± 0.5u 


3.1±0.6 C 


Q(14) 


0.6±0.3u 





Table A3. The comparison of our equivalent widths 
[mA] to data of van Dishoeck & de Zeeuw (1984) 
(D&Z) and Danks & Lambert (1983) (D&L) for 
HD 148184 of the (2,0) band. 



B(N") 


our results 


D&Z a 


D&L a 


P(6) 


0.7 ±0.4 


1.3 ±0.2 


0.96 


P(8) 


0.6 ± 0.4 




0.81 


P(4) 


1.5 ± 0.4 




0.44 


i?(10) 


0.7 ±0.4 




0.48 


R(2) 


2.7 ± 0.4 


1.7 ± 0.2 


1.33 


R(12) 




0.5 ±0.2 




R(0) 


1.4±0.4u 


1.4 ±0.1 


0.93 


Q(2) 


3.0 ±0.4 


2.4 ±0.1 


2.41 


fl(14) 




^ 0.5 




0(4) 


3.3 ±0.5 


2.5 ±0.2 


2.15 


P(2) 


0.5±0.5u 


^ 0.5 




Q(6) 


3.4 ±0.4 


2.4 ±0.2 


2.26 


i?(16) 




^ 0.5 




Q(8) 


1.7 ±0.4 


2.0 ±0.3 


1.67 


P(4) 


1.4 ±0.5 


0.8 ±0.2 


0.59 


Q(10) 


1.0 ±0.5 


0.8 ±0.2 


1.37 


P(6) 


1.0 ±0.4 


1.0 ±0.1 


0.67 


Q(12) 


0.6 ±0.4 


0.9 ±0.2 


0.78 


P(8) 




1.0 ±0.1 


0.89 


Q(14) 




0.6 ±0.1 


0.74 


P(10) 




0.5 ±0.2 


^ 0.67 


Q(16) 


0.8 ±0.4 







a Measurements were made with coude echelle spec- 
trograph (with Rcticon) which was fed by the 1.4-m 
coude auxiliary telescope at ESO, La Silla, Chile. 
The resolving power in those observations was 80,000. 



Table A5. The comparison of our equivalent widths 
[mA] to data of Gredel Sz Muench (1986) (G&M), van 
Dishoeck & Black (1989) (D&B) Crawford (1997) (Cr) 
and Gredel (1999) (G) for HD 169454 of the (2,0) band. 



B{N") 


our results 


G&M a 


D&B a 


Cr ft 


G c 


R(6) 


2.4 ±0.3 








< 2 


R(8) 


0.7 ±0.3 








1 


fl(4) 


3.7 ±0.3 


3.8 ±0.8 






4.1 ± 1 


R(10) 


0.4 ± 0.3u 










R(2) 


8.1 ±0.3 


6.1 ± 1.2 


6.0 ±0.5 




8.6 ± 1 


R(12) 


0.3 ±0.2 










R(0) 


8.5 ±0.3 


8.2 ±0.8 


5.6 ± 1.0 


6.23 ±0.50 


8.2 ± 1 


Q(2) 


10.3 ±0.3 


8.0 ± 1.6 


6.6 ±0.8 


8.14 ±1.18 


10.3 ± 1 


Q(4) 


7.9 ±0.3 


5.0 ± 1.0 


7.9 ±0.7 


3.98 ±1.03 


9.7 ± 1 


P(2) 


2.0 ±0.3 


2.3 ±0.5 


1.3 ±0.3 


2.82 ±0.55 


2.7 ±2 


Q(6) 


3.0 ±0.3 


3.1 ±0.6 


3.9 ±0.5 


2.89 ±0.59 


3.8 ± 1 


Q(8) 


2.3 ±0.3 


2.6 ± 1.0 


1.5 ±0.5 


1.28 ±0.59 


3.2 ± 1.5 


P(4) 


1.7±0.3 


2.0 ±0.8 


1.7 ±0.8 


1.68 ±0.62 


3.4 ± 1.5 


Q(10) 


0.8 ±0.3 


1.0 ±0.5 


1.0 ±0.8 




^ 1 


P(6) 


0.9 ±0.8 


1.9 ±0.8 


< 1.5 




3.0 ± 1 


Q(12) 


1.0 ±0.3 


1.0 ±0.5 


< 1.5 




^ 1 


P(8) 


0.8 ±0.3 










Q(14) 


0.6 ±0.4 










Q(16) 


0.6 ±0.4 











Measurements were made with coude spectrometer 
which was on the 2.1-m telescope at McDonald 
Observatory. The resolving power in those 
observations typically was 25,000 to 30,000. 
b P(4) and Q(8) lines arc blended (Federman et al. 1994). 
c Result is suspect because of possible stellar 
contamination (Fcderman ct al. 1994). 



a Measurements were made with coude echelle spectrograph (with R.cticon) 
which was fed by the 1.4-m coude auxiliary telescope at ESO, La Silla, 
Chile. The resolving power in those observations was 80,000. 
b The observations were obtained with the Ultra-High-Resolution Facility 



